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Abstract 


As  part  of  the  SYNOP  (Synoptic  Qcean  £rediction  experiment)  fidd  program,  twelve 
tall  moorings  measured  the  Gulf  Stream’s  temperature  and  velocity  fields  with  current 
metos  (CM)  at  nominal  depths  of  400  m,  700  m,  1000  m,  and  3500  m  f<»  two  years, 
from  May  1988  through  August  1990.  Simultaneously,  24  inverted  echo  sounders  (lES) 
monitored  the  thermocline  topography.  A  third  observational  component  the  ex¬ 
periment  was  the  release  of  isopycnal  RAFOS  floats;  70  such  floats  traversed  the  area 
monitored  by  the  CM  and  the  lES.  This  report  documents  the  methods  used  to  com¬ 
pute  vertical  motion  for  each  data  source,  and  the  differences  and  similarities  between 
the  three  methods.  Typical  vdocities  during  ‘strong’  events,  as  observed  by  or  inferred 
from  all  three  instruments,  was  1  —  2  mm  s~*  in  regions  near  the  center  of  the  Gulf 
Stream.  The  comparison  of  RAFOS  vertical  motions  and  vertical  motions  diagnosed 
from  CM  data  showed  excellent  agreement;  furtherm<»e,  CM  vertical  motions  and  lES 
vertical  motions  are  statistically  coherent  for  periods  longer  than  12  days.  We  conclude 
that  we  may  map  mesoscale  fields  of  t&(x,y,t);  the  fields  m<^>ped  are  consistent  with 
quasi-geostrophic  dynamics. 
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1  Introduction 


The  SXIioptic  (2cean  £re<lictioii  experiment  (SYNOP)  was  a  multi-institational  study 
of  the  Gulf  Stream  from  1987*1990.  D^gned  to  further  the  understanding  of  Gulf  Stream 
meanders,  and  to  facilitate  the  modeling  of  the  meanders,  the  SYNOP  experiment  was 
comprised  of  three  different  arrays:  an  Inlet  Array  near  Cape  Hatteras,  a  Central  Array 
near  68*W,  and  an  Eastern  Array  near  the  Grand  Banks  at  55*W.  Data  collected  in  the 
Central  Array  to  be  considered  in  this  report  came  from  12  tall,  high-performance  moorings 
on  which  were  four  current  meters  (CM)  (at  nominal  depths  of  400, 700, 1000,  and  3500  m), 
from  24  inverted  echo  sounders  (I^)  which  acoustically  monitored  the  depth  of  the  main 
thermocline,  and  from  approximately  70  isopycnal  RAFOS  floats  which  sampled  pressure 
and  temperature  along  an  isopycnal  surface  as  they  moved  downstream  throng  the  Gulf 
Stream. 

Vertical  motion,  to,  in  the  Gulf  Stream  has  been  observed  in  previous  studies  (e.g. 
Bower  and  RxMsby  1989)  and  inferred  in  others  (Hall  1986).  Here,  we  wish  to  compare  the 
vertical  motions  observed  by  the  isopycnal  floats  to  those  inferred  from  CM  obeecvaticms 
and  those  inferred  from  lES  observations.  The  three  methods  of  computing/determining 
vertical  motion  are  detailed  in  the  following  section,  with  23  case  study  events  following. 

j- 

2  Instruments  and  vertical  motion  calculations 

2.1  RAFOS  floats 

2.1.1  Methodology 

Bower  and  Rossby  (1989)  show  how  vertical  motion  is  related  to  RAFOS  float  motion, 
namely  that  upward  vertical  motion  occurs  as  the  float  moves  downstream  from  trough  to 
crest,  and  downward  motion  accompanies  motion  downstream  from  crest  to  trough.  The 
RAFOS  float  recorded  temperature,  pressure  and  location  at  eight-hour  intervals  as  it  moved 
through  the  stream.  The  method  for  determining  vertical  motion  from  these  variables  is 
straightforward.  The  computed  (3  ^)  is  estimated  using  a  second-order  centered 
finite  difference.  Because  St  =  S  hours,  the  vertical  motion  computed  is  an  average  over 
a  16-hour  interval.  Of  the  72  RAFOS  floats  that  were  released  as  part  of  the  SYNOP 
experiment  and  returned  useful  data,  more  than  90%  passed  directly  through  the  Central 
Array,  providing  a  large  dataset  of  observed  vertical  motion. 

2.1.2  Errors  in  Float  verticsd  motion  estimates 

Errors  introduced  by  the  float  not  being  perfectly  isopycnal  should  be  insignificant  for  the 
scales  of  motion  we  consider  (Rossby,  p«sonal  communication).  However,  because  floats 
move  through  features,  the  vertical  motion  measured  by  the  float,  a  sixteen-hour  average 
value,  may  differ  considerably  from  the  point  values  measured  by  CM  and  lES.  This  is 
most  obviously  the  case  if  the  float  passes  through  an  up/down  couplet  in  rixteen  hours,  in 
which  case  the  average  vertical  motion  may  be  near  zero,  although  the  values  of  upwelBag 
or  downwelling  within  the  couplet  (as  sampled  by  a  current  meter  mooring,  perhaps)  will 
not  be.  This  bias  error  may  underestimate  several  of  the  largest  w  events  by  as  much  as 
30%. 
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To  understand  this  bias,  consider  a  pattern  of  vertical  motion  that  along  a  fluid  particle 
trajectory  varies  sinusoidally  such  that  the  pressure  measured  by  a  RAFOS  float  and  the 
vertical  motion  will  be 

P  ^  A  sinojt  and  wtrue  =  A  u>  cosu/t,  (1) 

respectively.  How  is  the  computed  vertical  motion  affected  by  discrete,  rather  than  contin¬ 
uous,  sampling?  The  estimated  vertical  motion,  w,  is 

+  A01-i*n[w(t  -  At)]  stnwAt 

W 

The  ratio  approaches  unity  for  all  motions  that  have  period  T  much  longer  than  the 
sampling  time.  At,  but  it  crosses  zero  for  T  =  2At. 

The  RAFOS  floats  that  generated  the  data  used  in  this  study  sampled  data  at  At  =  8 
hours.  For  T  >  3  days,  for  example,  >  0.92;  however,  for  a  period  T  of  36  hours,  the 
estimated  vertical  motion  is  70%  of  true.  Vertical  motions  in  the  inertial  frequency  range 
(or  higher  frequencies)  are  significantly  underestimated  by  the  8-h  sampling  period.  For 
a  period  of  19  hours,  for  example,  the  estimated  vertical  motion  is  only  18%  of  the  true. 
However,  this  is  desirable,  because  we  are  using  the  floats  to  estimate  vertical  motions  on 
meso-  and  synoptic  scales. 

The  largest  vertical  motions  that  we  investigate  in  this  paper  are  typically  associated 
with  float  displacements  with  period  three  days  or  less.  Several  of  the  float  trajectories 
through  the  central  array  show  strong  up  and  down  couplets,  with  the  float  moving  np  and 
down  >  200  m  in  less  than  three  days.  Equation  (2)  indicates  that  the  vertical  motion  as 
measured  by  the  float  may  underestimate  the  true  vertical  motion  in  these  cases  by  around 
30%  if  the  dominant  period  of  the  float  is  about  36  hours. 


2.2  Current  meters 
2.2.1  Methods 

Bryden  (1976, 1980)  was  among  the  first*  to  use  the  backing  (turning  counterclockwise  with 
-'ecreasing  depth)  and  veering  (turning  clockwise  with  decreasing  depth)  of  current  meters 
to  infer  vertical  motion  in  the  ocean.  Following  his  notation,  and  that  of  Hall  (1986),  the 
temperature  equation. 


dT  ,  dT  dT  dB^  ^ 

”’■37=“' 

can  be  transformed,  using  the  thermal  wind  equation,  to 


'First  in  oceanography.  Arnason  (1942)  and  Panoflcsy  (1944)  developed  similar  methods  for  the 
atmosphere. 
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Here,  R  is  magnitude  of  the  current  velocity  (u^  +  po  is  mean  density  (1035  kg  m~^), 
g  is  gravitational  acceleration  (9.8  m  s~^),  /  is  the  Coriolis  parameter,  Og  is  the  effective 
thermal  expansion  coefficient  ^  +  |f  =  1.2  x  10“*A:<7  K~^)  appropriate  for 

the  T-S  relationship  in  mid-thermocline  in  the  Gulf  Stream  (Hall  1985),  ^  is  the  mean 
stratification  (0.02  K  m~^),  ^  is  the  angle  of  the  current  with  respect  to  east,  and  u  and  v 
are  the  eastward  and  northward  components  of  the  flow.  The  constants  used  in  this  work 
are  consistent  with  a  Brunt- Vaisala  period  of  about  twenty  minutes.  We  note  that 
could  be  biased  slightly  by  this  choice  of  constants. 

For  isopycnal  motions,  vertical  motion  is  proportional  to  the  two  terms  in  the  numerator 
of  (4)  or  (5).  The  first  describes  the  motion  of  the  isopycnals  themselves:  for  instance,  in 
the  absence  of  horizontal  motions,  if  the  temperature  at  a  current  meter  is  increasing, 
isopycnal  surfaces,  and  the  water  parcels  on  them,  must  be  descending.  The  second  term  in 
the  numerator  is  the  heat  advection  in  the  presence  of  sloping  isopycnal  surfaces,  which  is 
proportional  to  the  cross-frontal  component  of  velocity  in  the  presence  of  sloping  isopycnal 
surface,  as  in  (3),  or  equivalently  to  the  turning  of  the  current  vector  with  height,  as  in 
(5).  Where  currents  back  or  veer  with  height,  vertical  motion  must  exist,  if  the  temperature 
is  constant.  Consider  if  the  flow  at  400m  is  towards  the  southeast  at  80  cm  s~^and  the  flow 
at  1000m  is  towards  the  northeast,  also  at  80  cm  s~^.  For  this  case  the  thermal  wind  will 
point  due  south.  Cold  water  is  to  the  east,  and  warm  water  is  to  the  west,  of  this  vector, 
t.e.  isotherms  slope  up  to  the  east.  The  mean  eastward  flow  encountering  this  upward  slope 
to  the  isotherms  therefore  generates  upward  vertical  motion. 

2.2.2  Application 

The  CM  data  used  in  this  study  have  been  smoothed.  A  40-hour  low-pass  Alter  has  been 
applied  to  the  data  to  remove  high-frequency  signals  (diurnal  and  semi-diurnal  tides,  for 
example).  In  addition,  the  CM  data  have  been  adjusted  to  compensate  for  mooring  motion. 
When  strong  currents  surround  a  current  meter  mooring,  it  will  tilt  over  more  than  when 
in  an  environment  of  weaker  currents.  Sensors  therefore  do  not  always  remain  at  the  same 
level.  The  method  of  Hogg  (1991)  interpolates  (or  extrapolates)  the  current  meter  data 
at  variable  levels  to  constant  horizons,  in  this  case  at  400,  700,  and  1000  m  below  the  sea 
surface.  We  computed  the  vertical  motion  using  data  with  and  without  the  mooring  motion 
correction  and  noted  negligible  differences.  Because  the  corrected  data  also  yield  vertical 
motions  at  the  same  level  throughout  the  domain  we  used  in  this  study  data  corrected  for 
mooring  motion. 

The  vertical  motion  was  computed  using  (4),  with  the  vertical  derivatives  (^  and 
In)  approximated  as  centered  finite  differences  using  data  at  400  and  1000  m.  Current 
vdodties  (u  and  v)  were  averages  of  values  at  400  and  1000  m.  Temperature  tendencies, 
(«).  were  computed  using  a  second-order  centered-in-time  finite  different?  scheme,  with 
At  =  12  hours,  and  data  at  700  m.  The  vertical  motion  computed  can  be  considered  as  a 
representative  value  about  halfway  between  the  400  and  1000  m.  We  have  used  the  formula 
including  -  «|n  for  two  reasons.  The  data  includes  u  and  v  at  three  levels,  so  teas  data 
processing  is  necessary  than  if  we  used  the  formula  including  ^  (although  interpretation 
of  the  results  is  easier  using  the  turning  current  vector  with  hei^t).  Furthermore,  the 
term  introduced  spikes  into  the  computed  vertical  motions  when  current  speed  profile  was 
non-linear. 

Because  the  current  speeds  used  in  (4)  are  at  two  levels,  400  and  1000m,  we  originally 
computed  the  temperature  tendency  term,  using  the  average  of  the  values  at  1000  and 
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400in  to  estimate  the  water  column  average.  This  introduced  an  error,  however,  for  those 
cases  in  which  the  400-m  CM  was  above  the  main  thermocline  and  the  1000-m  CM  was 
bdow.  Both  current  meters  then  recorded  small  temperature  tendencies,  even  though  the 
movement  of  the  main  thermocline  was  large  (as  measured  by  the  temperature  change  at 
700  m).  For  an  example,  see  Case  1,  section  3.1.  Although  using  the  averaged  temperature 
tendency  caused  an  error  as  just  described,  no  similar  error  is  evident  from  using  the  average 
of  tt  and  V  values  at  400  and  1000  m  to  represent  the  flow  at  700  m. 

2.3  lESs 

Inverted  echo  sounders  sit  on  the  ocean  floor  and  acoustically  monitor  the  depth  of  the  main 
thermocline,  denoted  as  the  depth  of  the  12^  isotherm,  i.e.  zij  (Watts  and  Rossby  1977). 
Vertical  motion  is  computed  from  lESs  by  using  objectively  mapped  lES  data  (Tracey  and 
Watts  1991)  to  determine  a  streamfunction  and  therefore  a  vorticity  field  as  in  Kim  (1991). 
As  with  CM  data,  all  lES  data  are  40-hour  low-pass  filtered  to  remove  high-frequency 
signals.  This  smoothing  is  done  before  the  objective  mapping.  The  lES  i)dds  are  objectively 
mapped  daily  for  26  months,  yielding  estimates  not  only  of  streamfunction  and  vortidty, 
but  of  the  time  tendencies  of  both,  if  the  motion  is  assumed  to  be  quasi-geostiophic,  then 
the  vorticity  equation  as  in,  for  example,  Holton  (1979), 

^  •  ViCg  +  /)  +  (*) 

can  be  rearranged  to  determine  vertical  stretching.  Vertical  motion  in  the  Gulf  Stream 
has  primarily  a  first  barodinic  mode  structure  (Hall  1986;  Rossby  1987).  The  maximum 
vertical  motion  occurs  below  the  main  thermo^ne  (Gill  1982,  fig.  6.14);  even  for  y«y 
shallow  msdn  thermodlnes  (t.e.,  Zi2  around  200  m  -  much  shallower  than  for  any  cd  the 
23  cases  considered  here),  the  methods  used  by  Piclcart  and  Watts  (1990)  show  that  the 
maximum  in  vertical  motion  is  below  500  m.  Vertical  motion  in  the  the  upper  ocean  at 
depth  z  can  be  approximated  as 

^itts  =  ’  V(Cp  +  /)]  (7) 

using  data  at  400  m. 

We  have  used  this  method  of  computing  vertical  motions  instead  of  the  quasi-geostrophic 
u;-equation  as  used  by  Leach  (1987),  Tintor4  et  al.  (1991)  and  Pollard  and  R^ier  (1992) 
because  it  is  more  suited  to  our  data.  In  the  other  studies,  roughly  83moptk  CTD  data  were 
available  to  provide  data  at  different  levels  as  is  necessary  for  Q-vector  computatkma.  We 
have  a  streamfunction  map  at  just  one  level,  but  at  a  sequence  of  times,  making  (6)  a  more 
natural  way  to  compote  vertical  motion.  Our  assumption  of  primarily  first-barodiuic  mode 
structure  is  consistent  with  the  definition  of  the  geostrophic  streamfunction  (Kim  1991), 
namely 


Vg=:^kx  v*ia  ,  (8) 

Jo 

where  we  have  used  for  p*(s  g^)  a  value  that  is  a  weak  function  of  z\2.  The  geostrophk 
barodinic  current  from  which  we  compote  ^  is  valid  at  400m,  and  is  referenced  to  3500m, 
where  no  motion  is  assumed. 
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To  generate  the  field  of  vertical  motion,  the  time  tendencji  ^|tsto  estimated  using 
a  centered  finite  difference  ( with  At  =  24  hours.  Vorticity  advection 
was  computed  at  time  to  with  horizontal  derivatives  computed  using  second-order  centered 
finite-differencing.  We  have  smoothed  the  lES  streamfunction  fields  with  three  applications 
of  a  Shapiro  second-order  filter  (Shapiro  1970)  to  reduce  the  amount  of  noise  in  the  highly 
differentiated  vertical  motion  field.  Such  smoothing  reduces  values  of  extrema  in  the  resul¬ 
tant  vertical  motion  field  without  changing  their  location.  The  filter  wavelength  cutoff  is 
—6  dB  at  133  km;  features  with  lengthscales  =  k~^  >  37  km  are  passed  with  greater  than 
80%  smiplitude.  The  gridded  values  of  interpolated  to  each  CM  location  so  that 

direct  comparisons  could  be  made  between  and  Wjgg. 

The  vertical  motion  as  diagnosed  by  (7)  is  forced  by  two  effects,  the  local  change  and 
the  advective  change  and  Vg  •  v(Cs  +  /),  respectively).  Advection  of  vorticity  causes 
vertical  motion  because,  as  explained  in  Holton  (1979),  the  vertical  motion  is  required  to 
keep  the  ocean  hydrostatic  and  geostrophic.  Specifically,  vertical  motion  acts  to  depress 
the  thermocline  (for  the  case  of  anticyclonic  vorticity  advection  in  which  case  a  crest  is 
approaching)  or  to  raise  the  thermocline  (for  the  case  of  cyclonic  vorticity  advection,  asso¬ 
ciated  with  an  approaching  trough).  If  the  vorticity  locally  is  changing  in  a  positive  sense, 
t.e.,  becoming  more  cyclonic  with  time,  then  we  should  expect  downward  motion:  vortidty 
is  most  easily  produced  locally  by  convergence,  which  will  be  accompanied  by  downward 
motion.  Note  that  if  the  advection  at  a  point  is  cyclonic  and  the  local  tendency  is  positive 
(i.e.,  becoming  more  cyclonic  with  time),  the  two  effects  compete.  If  the  local  tendency  is 
exactly  equal  to  the  advective  tendency,  then  there  is  no  vertical  motion  (and,  inddentaHy, 
the  phase  speed  and  the  float  parcel  speed  are  equal).  Normally,  the  advective  tenn  b 
larger;  that  is  the  case  for  75%  of  the  23  cases  considered  here. 

Finally,  note  that  the  mapped  lES  fields  are  a  refined  version  of  those  described  in 
Tracey  and  Watts  (1991);  in  those  fields,  the  first  guess  field  (t.e.,  the  so-called  ‘^ean’* 
field)  was  a  broad  diffuse  Gulf  Stream  (see  their  Figure  6).  The  fields  in  Tracey  and  Watts 
(1991)  contained  measurable  bias  that  could  be  traced  to  this  first-guess  mean  field.  The 
first-guess  field  used  for  this  study  is  a  31-day  running  mean  field  produced  from  the  data  in 
Tracey  and  Watts.  Because  the  first-guess  field  is  closer  to  the  actual  field,  the  objectively 
analyzed  fields  are  more  accurate.  A  further  refinement,  that  used  pseudo-IES  data  from 
CM  data  to  fill  in  gaps  in  the  lES  data,  and  for  which  lES  z\2  data  was  recomputed  to  reflect 
new  B-intercepts  based  on  a  different  Q  curve  (see  Watts  and  Rossby  1977  for  definitkms), 
has  not  been  incorporated  into  this  study. 

3  Case  Study  Comparisons 

To  compare  with  we  selected  those  floats  that  passed  within  10  km  of  a  mooring 
on  which  all  of  the  top  three  CMs  were  functioning.  The  10-l;m  cutoff  was  chosoi  as 
a  compromise  between  increasing  error  (due  to  the  highly  sheared  environment)  as  the 
distance  increased  from  CM  to  float  versus  obtaining  a  useful  number  (23)  of  comparisons. 
The  requirement  that  the  three  top  CM  functioned  was  to  minimize  error  in  thb  prdiminaiy 
study:  as  previously  noted,  the  vertical  motion  could  be  calculated  using  data  at  only  2  CM 
levels,  but  this  introduces  an  error;  furthermore,  requiring  the  presence  of  three  bVds  of 
data  minimizes  the  error  in  the  mooring  motion  correction  algorithm.  Note  that  the  data 
are  chosen  based  on  the  CM;  these  sites/times  chosen  may  not  have  optimal  lES  coverage. 
We  will  note  where  this  may  cause  errors. 
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We  carefully  scrutinized  the  CM  data  records  to  make  certain  that  submesoscale  coherent 
vortices  (SCVs),  as  described  in  Bane  tt  al.  (1989),  were  not  present  when  comparisons 
were  ma^.  Because  SCVs  typically  affect  only  one  CM  on  a  mooring,  canung  up  to  a 
360*  rotatkm  in  the  current  vector  as  wdl  as  a  temperature  perturbation,  their  presence 
obviously  would  adversely  impact  this  primarily  quasi-geostrophic  study.  No  evidmioe  ci 
obvious  SCVs  was  found  in  any  5-day  interval  centered  on  the  23  closest-passing  float  cases 
studied  here.  It  is  possible,  however,  that  there  could  be  cases  in  which  subtle  submesoscale 
features  partially  detected  in  CM  data  and  not  easily  identifiable  as  SCVs  (and  therdbre 
present  in  tOcw)  filtered  out  of  the  objectively  mapped  lES  streamfunction  Adds  (and 
therefore  not  present  in  iC/cs),  t.e.  subm^oscale  eddies  could  cause, errors  as  shown  in 
Panofsky(1951). 

For  each  of  the  case  studies  that  follow,  several  figures  are  presented.  The  track  of  the 
RAFOS  float  through  the  Central  Array,  as  well  as  the  pressure  record  of  the  float  will 
be  shown.  Mapped  fields  of  10, will  also  be  shown,  as  will  mapped  fields  of  ^  and 
^  '  V(C«  +  /)•  >•«•,  the  two  components  of  We  will  also  show  a  figure  the  torms 
contributing  to  Wcu- 

3.1  Case  1:  Float  123,  CM  H4,  4  July  1088,  0000  UTC 

3.1.1  Overview 

The  time  surrounding  float  123*^  passage  by  CM  H4  (Figure  1)  is  one  during  whidi  a  series 
of  crests  and  troughs  propagated  steadily  through  the  Central  Array.  This  Is  leflectsd  fai 
the  float’s  motion  through  several  crests  and  troughs  (made  obvious  by  the  large  vertical  ex¬ 
cursions  associated  with  changes  in  path  curvature),  in  the  IBS  *12  topop»phy  (FIgare  2a), 
which  shows  a  crest  (at  s  ==  -100)  and  a  trough  (at  z  =  40)  in  the  CoitrM  Array,  ai^ 
in  the  CM  records,  which  show  a  distinct  oscillation  in  current  direction  and  temperature 
(Figure  3a,b)  caused  by  propagating  features. 

3.1.2  Vertical  motion  results 

=  -0.58  mm  s”*  ,  Wjgg  =  -1.42  mm  s“*  ,  =  -0.98  mm  s"* 

Float  123  moved  from  532  to  550  to  565  db  (Figure  1)  for  the  three  times  closest  to  its 
pass  by  CM  H4,  a  mean  vertical  motion  of  -0.58  mm  s~*  for  those  16  hours.  When  float 
123  was  closest  to  CM  H4,  the  CM  data  suggests  that  vi^rous  downwelflng  was  occuniiC. 
Temperatures  were  increasing  as  isotherms  descended  (Figure  3b)  in  advance  at  the  ap- 
proaching  crest  (at  —1.13  mm  s'*  )  concomitant  with  a  slight  veering  in  the  currents  with 
decreasing  depth  (Figure  3a)  (suggesting  weak  warm  advection  and  upward  vertical  motion 
at  0.15  s~*  (Figure  3e)).  Note  that  for  this  case,  using  the  average  at  the  temperature 

change  at  400  and  1000m  to  estimate  the  temperature  change  at  700m  yMds  an  incorrect 
answer  (Figure  3d4)-  Wjgg  for  the  same  time  (0000  UTC  4  July  1988)  was  — 1 JIO  nun  s***  . 
Based  on  adve^ion  alone  (Figure  2d),  vertical  motion  should  be  downward  at  1.67  nun  s~* ; 
strong  anticydonic  vortidty  advection  forces  the  local  vortidty  to  beoNne  morn  and  mote 
nq^ative.  Anticydonic  spin-ap  would  be  accompanied  with  downward  motion  to  force  down 
the  thermocline.  The  local  tendency  at  the  H4,  however,  amtributed  to  upward  vertical 
motion  (0.25  mm  s"*  ).  Upward  motion  is  forced  by  divergence,  which  destroys  vortidty. 
This  is  an  example  of  the  advective  part  of  the  vortidty  equation  being  counteracted  by 
the  local  tendmicy. 
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Lonfitud*  (W) 

Figure  1:  Preerare  record  in  db  (top)  and  borixontai  track  (bottom)  -  poiata  after  tbe  iiet 
(lavnch)  point  are  ei|^t  boon  ^>art  -  cS  EAFOS  float  123.  Clf  locatkma  wHUn  tbe  Central 
Array  are  indicated  by  open  drclee.  Tbe  three  doeest  RAFOS  fixes  to  CM  H4  are  bonad. 
Lanncb  date  of  float:  1037  UTC,  29  Jane  1988. 
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Figure  2:  (A)  Objectivdy  mapped  depth  of  12*  isotherm  (sia)  as  measured  Iqr  lES  data  for 
0000  UTC  4  July  1988,  the  time  when  float  123  was  closest  to  CM  H4.  Cmitoar  interval 
every  100m.  AIm,  deq>  current  stick  vect<ws  plotted  at  CM  locatkms,  with  a  hqr  (vahw 
as  0.25  m  s~')  plotted  at  x  s  -155,  y  s  -IM.  (B)  Objectivdy  mapped  xia  topopaphy 
(Ixdd  solid  lines,  oontonrs  at  150,  350,  550,  750,  and  950  m)  and  portion  of  D  farced  by 
vortidty  advecthm  (thin  sdid  lines,  contonr  interval  0.50  mm  s'* ,  negative  valnes  dished). 
(C)  Objectively  mi^iped  xia  topography  (bold  solid  lines,  conUmrs  at  150,  350,  550,  750, 
and  95te)  and  portioa  of  D)  forced  by  local  vwticity  tendenqr  (thin  lines,  contonr  i^srvnl 
0.25  mm  s'** ,  negative  valnes  dadied).  (D)  Objectii^y  mapp^  xia  topogr^hy  (bdd  sold 
fines,  oontonrs  at  150, 350,  550,  750,  and  950m)  and  total  vertical  motion  competed  from 
(7)  n^ig  lES  data  (thin  fines,  contonr  interval  0.5  mm  s~*  ,  negative  valnes  dashed).  Ihr 
all  four  plots,  the  trade  of  RAFOS  float  123  is  indicated  by  dots,  and  the  location  vt  CM 
H4  is  marked  Inr  a  drde. 
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Figure  3:  CM  H4  daU  for  5  dayi  surrounding  float  123’s  closest  approach  to  H4,  wkiA  time 
is  indicated  by  a  vertical  line.  (A)  Direction  (in  *  )  towards  wbi^  tlM  cuieat  is  ioerfaig 
(sdid  line:  400m;  kmg  dashed  Une:  700m;  short  dashed  line:  1000m).  (B)  Thmpentwe 
at  CM  site  in  degrees  cdsius  (sdid  Une:  400m;  bng  dashed  line:  700m;  sbort  dashed  Bam 
1000m).  (C)  Current  Speed  in  m  8~*(solid  line:  400m;  long  dashed  line:  TOOnt;  short  dashed 
line:  1000m).  (D)  Vertical  motion  (in  mm  s~*  )  finrced  by  local  tanperature  tsndsnqr  (soBd 
line:  using  the  observed  tonpoature  toidency  at  700m;  dashed  line:  udag  the  aeeracs  of 
400  and  1000  m  temperature  tendencies  to  estimate  the  700m  tonperataie  tsndsnqr).  (E) 
Vertical  motion  (in  mm  s~*  )  forced  by  baddng/veeriag  current  (dashed  Base  ooaiv*^ 
as  in  (4)  n^g  -  «|||;  solid  line:  computed  as  in  (5)  using  (F)  IMal  vertical 

motkm  in  mm  ,  t.e.  the  sum  ci  E  (dashed  line)  and  E  (solid  and  dashed  Baa  as  in  D). 
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3.2  Case  2:  Float  120,  CM  14,  23  September  1088,  0800  -  1000  UTC 
3.2.1  Overview 


RAFOS  float  120  spent  considerable  time  in  the  Central  Array  (Figore  4).  Althoogh  it 
entered  moving  steadily  on  an  anticyclonically  curved  path,  it  then  ^>paiently  fc^p*** 
temp<mtfily  from  the  stream  and  moved  only  slowly  before  a  meander  trouf^  moved  thrcragh 
and  started  to  steer  the  float  again  steadily  to  the  northeast.  The  lES  U^tography  for  the 
dates  when  the  float  was  closest  to  CM  14  (Ilgnres  5, 6)  do  show  a  considerable  trongh  near 
the  center  cS  the  Central  Array.  CM  curraits  (Figure  7)  ate  weak  (less  than  0.5  m  s'*). 
Temperatures  are  slowly  increasing. 


Figure  4:  As  in  Figure  1,  but  for  float  129.  The  four  closest  RAFOS  flam  to  CM  14  ate 
boomd.  Launch  date:  0838  UTC,  31  August  1988. 


3.2.2  Vertical  motion  results 

w  0.19  mm  s'*  ,  »  0.53  mm  s'*  ,  w,gg  s  0.04  mm  s'*  (0800  UTC) 
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Figure  5:  Ai  in  Figure  2,  except  for  0800  UTC  23  September  1988.  The  track  of  font  129 
is  indicated  by  dots,  and  the  location  of  CM  14  is  marked  by  a  circle. 


^HAr  —  ®  •"*  (1600  UTC) 

Fev  the  four  dosert  float  locations  to  CM  14  (Figure  4),  the  float  npweOed  fitom  717  to 
711  to  706  to  669  db  (between  00  UTC  23  Sep  and  00  UTC  24  Sep);  the  two  wmam.  rahMS 
of  vwtical  mothm  are  0.19  and  0.72  mm  s~*  .  Clearly,  the  float  is  b^inning  to  acedsrate 
upwards  as  it  moves  out  its  position  at  the  base  of  a  trough*  CM  data  also  suggest  risfaig 
motion.  Temperatures  are  coding,  suggesting  that  isotherms  are  rising  (Figure  7d^).  At 
the  same  time,  there  is  sigiuficant  veering  to  the  currents  from  whidi  we  caa  iufcr  iWug 
mdkm  (Ftgure  7a-c).  The  vertical  motkm  as  estimated  by  the  thermocBne  topogn^  ^  (lES 
data)  is  consideraMy  smaller,  about  0.04  mm  s~*  for  both  times  (Figures  5,6).  Thai  pwt 
due  to  adveetkm  is  nearly  sero  and  the  part  due  to  local  tendency  is  about  0.04  mm  s~*  . 
Ws  note,  however,  that  a  large  area  at  ascrat  is  located  just  to  the  mwth  at  14. 


wm  (Undmcy)  wng  (total) 

flgnra  0:  Ai  ia  Flgua  5,  except  for  1600  UTC  23  September  1988. 
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I  Figure  7:  As  in  Figure  3,  except  for  CM  14  for  the  live  days  surrounding  float  129*s  doeeet 

^proachee  to  14,  both  of  which  times  are  indicated  by  vwtical  lines. 

I 
I 
I 
I 
I 
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3.3  Caae  3:  Float  138,  CM  14,  8  December  1988,  1600  UTC 
3.3.1  Overview 

The  Golf  Stream  in  the  period  daring  which  float  136  (Figare  8)  moved  through  the  Central 
Array  was  characterised  by  a  succession  of  crests  and  meanders.  Float  136  initially  moved 
fairly  slowly  through  the  Central  Array  along  the  southern  edge  of  the  stream.  Halfway 
through  the  array,  however,  the  float  accelerated  to  the  northeast,  M>pvently  influoiced  by 
an  approaching  propagating  trough  (Figare  9). 


Figare  8:  As  in  Figure  1,  but  for  float  136.  The  three  closest  RAFOS  Axes  to  CM  14  are 
boxed.  Launch  date:  0921  UTC,  19  November  1988. 


3.3.3  Vertical  motion  results 

=  0.81  mm  s“*  ,  =  0.43  mm  s”*  ,  Wjg,  =  0.17  mm  s"* 

Float  136  moved  from  669  to  668  to  622  db  during  the  time  it  was  closest  to  CM  14  (Fig¬ 
ure  8),  reflecting  the  acceleration  up  and  northeastward  out  of  a  meandpr  trough.  CM  data 
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Figure  9:  As  in  Figure  2,  except  for  1600  UTC  8  December  1988.  The  track  ct  float  136  b 
indicated  by  dots,  and  the  location  of  CM  14  is  marked  by  a  circle. 


at  the  same  time  show  currents  that  are  veoing  with  decreasing  depth,  suggesting  upward 
motion  (Figure  10a,c,e).  At  the  same  time,  CM  temperatures  are  steady  (Figure  10b);  the 
majority  of  vwtical  motion  diagnosed  from  CM  data  is  related  to  the  change  of  the  current 
vectw  with  height.  Vertical  motion  determined  from  TES  data  shows  weak  upward  motion 
forced  by  vortidty  advection  (Figure  9b),  and  weaker  still  downward  motion  rdated  to  the 
vortidty  tendency  (Figure  9c);  the  net  motion  (Figure  lOd)  is  upward. 
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3.4  Caae  4:  Float  141,  CM  14,  2  November  1988,  1600  UTC 
8.4.1  Overview 

RAFOS  float  141  (Figure  11)  moved  steadily  through  the  Central  Array  at  a  time  during 
which  flow  was  from  the  northwest  to  the  southeast,  suggesting  a  large  ridge  upstream  of  the 
Central  Array.  However,  there  are  wiggles  in  the  float  path  suggesting  smaller-scale  crests 
and  troughs  are  propagating  through  the  iarger-scale  ridge,  and  this  view  is  confirmed  by 
the  lES  Zi2  topography  (Figure  12). 


n&Yi  ama  uuMCB 


Longitude  (W) 

Figure  11:  As  in  Figure  1,  but  for  float  141.  The  three  closest  BAFOS  fixes  to  CM  14  are 
boxed.  Launch  date:  0921  UTC,  28  October  1988. 


3.4.2  Vertical  motion  results 

=  -0.87  mm  s“*  ,  =  -0.94  mm  s~*  ,  w,gg  =  -0.43  mm  s“* 

RAFOS  float  141  downweUed  from  589  to  603  to  639  db  (Figure  11)  in  the  period  closest 
to  CM  14,  which  is  a  mean  motion  of  -.87  mm  s"^  .  CM  data  for  the  time  of  closest  ^>- 
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Figure  12:  As  in  Figure  2,  except  for  1600  UTC  2  November  1988.  The  track  ct  float  141 
is  indicated  by  dots,  and  the  location  of  CM  14  is  marked  by  a  circle. 

proach  also  suggests  downwelling;  although  the  currents  are  veering,  suggeatlag  npuralUag 
(Figures  13a, c,e),  temperatures  are  wanning  rapidly  enough  that  the  accompanying  down- 
welling  dominates  (Figures  13b,d4')-  Note  that  here,  again,  using  the  average  oi  400  and 
1000  tmnperatore  tendencies  to  predict  the  tendency  at  700  m  leads  to  an  underpredictkMi 
of  the  magnitude  of  the  vertical  motion  caused  by  the  local  temperature  tendency,  and 
hence  a  sign  error  in  the  resultant  total  vertical  motion  (Figures  13d,f).  lES  data  shows 
that  CM  14  is  downstream  of  a  crest;  anticyclonic  vorticity  advectkm  associated  with  this 
crest  is  forcing  downward  motion  (Figure  12). 
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Figure  13:  As  in  Figure  3,  except  for  CM  14  for  the  five  days  surrounding  float  Idl’s  doeest 
approaches  to  14,  which  time  is  indicated  by  vertical  line. 
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3.5  Case  5:  Float  175,  CM  I4,  20  January  1089,  1600  UTC 
S.5.1  Ovarview 

As  was  the  case  for  float  141,  EAFOS  float  175  traversed  the  Central  Array  moving  from 
northeast  to  southwest  on  a  track  consistent  with  an  upstream  crest  (figure  14).  There 
were  wiggles  in  the  track  that  suggest  smaller  scale  features  propagating  through  the  larger 
scale  crest.  This  is  confirmed  by  the  mr^)  of  xu  (Figure  15). 


figure  14:  As  in  Figure  1,  but  for  float  175.  The  three  closest  RAFOS  fixes  to  CM  14  are 
boxed.  Launch  date:  1005  UTC,  14  January  1989. 


3.5.2  Vertical  motion  results 

=  -0.56  mm  8“*  ,  =  -0.05  mm  8~*  ,  XD,gg  *  -0.63  mm  s"* 

RAFOS  float  175  moved  in  the  vertical  from  328  to  364  to  361  db  during  the  time  when 
it  was  closest  to  CM  14  (Figure  14),  a  mean  downwdling  of  —0.56  mm  s~*  .  Note  that 
float  is  near  the  bottom  of  a  trough  and  is  about  to  commence  upwdling  downstream  of  it. 
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Figure  15:  As  in  Figure  2,  except  for  1600  UTC  21  January  1989.  The  track  of  float  175  is 
indicated  by  dots,  and  the  location  of  CM  14  is  marked  by  a  circle. 


CM  vertical  motion  is  also  flipping  sign  (Figure  16)  as  the  float  passes  by  14,  reflecting  the 
propagation  of  the  meander  crest.  Vertical  motion  is  stiU  slightly  downward,  however,  as 
measured  by  the  CM  when  the  float  is  close  by.  Vertical  motkm  as  diagnosed  by  vertical 
stretching  data  at  400  m  is  also  downward;  CM  I4  is  downstream  of  a  crest  (Figure  15). 
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Figun  16:  At  in  Figon  3,  except  for  CM  14  for  the  five  days  tvrrouuliag  float  175*t  doteet 
approechet  to  14,  which  time  it  indicated  by  vertical  line. 


3.6  Caa*  6:  Float  176,  CM  IS,  7  Fabmary  1669, 0000  UTC  -  8  Fabraary  1669, 
0000  UTC 

3.6.1  Ovwviaw 

Float  176  inovad  alowly  through  the  Ceatral  Array  at  first  (I%are  17),  tradag  oat  aa 
aatiqrdonic  path,  b^ore  being  expelled  from  the  stream,  after  whkh  time  it  drifted  very 
slowly  southward.  lES  data  (Figores  18-  21)  show  the  crest  on  the  eastern  edge  of  the 
Central  Array,  and  also  suggests  a  sharp  southern  edge  to  the  stream:  float  176  is  moving 
very  slowly  in  a  region  that  is  dose  to  the  main  part  ot  the  Gulf  Stream. 


Longitude  (W) 


Figure  17:  As  in  Figure  1,  but  for  float  176.  The  seven  doeest  RAFOS  fixes  to  CM  15  are 
boxed.  Launch  date:  1131  UTC,  25  January  1989. 


3.6.2  Verticnl  motion  resuUa 

s  -0.04  mm  -  0.02  mm  s'* ,  w,gg  =  0.25  mm  s~*  (0000  UTC  7  FUnuary) 
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Figare  18:  As  in  Figure  2,  except  for  0800  UTC  7  Febmery  1989.  The  tredt  of  ioet  176  is 
indketed  by  dots,  end  the  bcatioa  of  CM  15  is  merked  by  »  circle. 


s  -0.05  mm  s'** ,  *  ^iss  ~  (0800  UTC  7  FWmmry) 

a  0.04  mm  s'* ,  m  0.02  mm  s'* ,  w,„  a  0.27  mm  s'*  (1600  UTC  7  Ibbruiy) 
a  -0.03  mm  s'*  ,  a  -0.05  mm  s'*  ,  w,„  a  0.26  mm  s'*  (0000  UTC  8  fMin> 

•*y) 

As  the  RAFOS  Hoot  drifts  soathward  ia  the  vidnity  of  CM  15  (H  leoMias  witkia  10  km 
15  for  a  day),  it  stays  betweea  761  aad  766  db  (Figure  17),  so  vertical  awtioa  is  as^gfola.  fti 
additkw,  CM  data  disgaoees  small  vertical  motkm  (Figure  22);  there  Is  veariag  aad  baddag 
<A  cnrreats,  but  the  cuneats  are  extroady  weak.  Similady,  local  temperature  dMuigm  are 
very  saialL  lES  data  (Figures  18-  21),  however,  diagaoses  upward  vertical  aiotioB  forced  by 
cydoak  vortidty  advectimi  (local  teadeacy  tA  vortidty  is  aear  sero).  is  aeady  0.35 
aim  s'*  for  tbe  dayloag  period  whea  the  float  is  close  to  15;  however,  15  is  aear  tbe  edge  of 
the  lES  aiappiag  rsgioa,  aad  the  derivatives  required  to  compute  verity  aad  its  advectioa 
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Ilgan  19:  As  ia  flgara  18,  except  kt  1600  UTC  7  Febnury  1989. 
may  atrt  be  accurate. 
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Figu*  20:  At  in  flgun  18,  •xcq>t  for  0000  UTC  8  February  1989. 


I  flgon  21:  As  in  Figore  18,  except  for  0800  UTC  8  Fsbraarj  1089. 

I 

I 

I 
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Figure  22:  As  in  Figure  3,  except  for  CM  15  for  the  five  days  sorroanding  float  176*s  dosest 
^»proaclMB  to  15,  which  times  are  indicated  by  vertical  lines. 


3.7  Case  7:  Float  194,  CM  IS,  IS  November  1989,  0000  UTC 
3.7.1  Overview 

BAFOS  float  194  moved  through  the  central  array  during  a  time  of  ngnificant  stream¬ 
ring  interaction.  This  is  reflected  in  its  unusual  south-to-north  path  through  the  array 
(Figure  23),  moving  around  a  cyclonic  add  eddy  in  the  southern  part  oi  the  domain,  aa<* 
around  a  crest  in  the  northern  part  of  the  domain  (Figure  24). 


I  ■  1  1  ■  I  •  1.  I  ■  1  -  I  .  I 

77  76  75  74  73  72  71  70  69  68  67  66  65 
Longitude  (W) 


Figure  23:  As  in  Figure  1,  but  for  float  194.  The  three  closest  BAFOS  fbms  to  CM  IS  an 
boxed.  Launch  date:  1843  UTC,  5  November  1989. 


3.7.2  Vertical  motion  results 

~  *  “'cv  “  »“*  *  “>/»s  =  -0.22  mm  s"* 

As  float  194  moved  around  the  odd  eddy  past  15,  it  upwdled  from  663  to  662  to  622 
db  (Figure  23),  a  mean  value  of  .71  mm  s~*  .  This  is  just  as  expected  from  the  work  of 
Bower(1989);  the  float  has  emerged  from  a  trough  and  is  moving  towards  a  crest.  The  path 


Figure  24:  As  in  Figure  2,  except  for  0000  UTC  15  November  1989.  The  track  of  float  194 
is  indicated  by  dots,  and  the  location  of  CM  15  is  marked  by  a  circle. 


is  more  convoluted  than  normal,  however.  CM  data  at  the  same  time  at  15  also  suggest 
upwdling  (Figure  25):  Currents  are  veering,  and  temperatures  are  cooling  slightly.  Both 
of  these  win  force  upward  vertical  motion.  lES  data  (Figure  24)  for  this  time,  however, 
suggests  downward  motion.  Both  vortidty  advection  and  local  tendency  forcing  down¬ 
ward  motion.  A  passible  cause  for  the  discrepancy  in  diagnosed  vertical  motion  is  discussed 
in  the  next  section. 
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Figure  25:  As  in  Figure  3,  except  for  CM  15  for  the  five  days  surrounding  float  lfl4*s  closest 
approaches  to  15,  which  time  is  indicated  by  vertical  line. 


3.8  Caie  8:  Float  194,  CM  14,  16  November  1980,  0000  UTC 
3.8.1  Overview 

RAFOS  float  194’8  unasual  path  (Figure  26)  south  to  north  from  ring  to  crest  (Figure  27) 
is  described  in  Section  3.7. 


nmAmuimai 


Longitude  (W) 

Figure  26:  As  in  Figure  23.  The  three  closest  RAFOS  fixes  to  CM  14  are  boxed.  Launch 
date:  1843  UTC,  5  November  1989. 


3.8.2  Vertical  motion  results 

^RAr  —  0'®®  »  ^cM  —  ^*20  mm  s'*  ,  Wjgg  =  -1.06  mm  s“* 

Float  194  continued  upwelling  as  it  passed  CM  14,  moving  from  619  to  584  to  570  db 
(Figure  26).  This  is  a  mean  upwelling  of  0.86  mm  s~*  .  CM  data  at  14  also  diagnosed 
upward  motion  (Figure  28),  forced  mostly  by  strong  cooling  observed  at  the  current  meter 
site  (Figure  28b, d).  Vertical  motion  caused  by  currents  veering  with  bright  is  smaller  in 
size.  In  contrast  to  the  upwelling  measured  by  the  float  and  diagnosed  by  the  CM,  I£S  data 
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Figure  27:  As  in  Figure  24,  except  for  0000  UTC  16  November  1989.  The  track  of  float  194 
is  indicated  by  dots,  and  the  location  of  CM  14  is  marked  by  a  circle. 


(Figure  27)  suggests  downwelling.  Weak  downwelling  is  being  forced  by  vortidty  advectioa, 
supplementing  the  stronger  downwelling  forced  by  the  local  change  in  the  vortidty.  This 
considerable  disagreement  between  W/gg  and  the  other  two  methods  is  discussed  in  the  next 
section. 
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Figure  28:  As  in  Figure  3,  except  for  CM  14  for  the  five  days  surrounding  float  194*s  closest 
approaches  to  14,  which  time  is  indicated  by  vertical  line. 
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S.9  Case  0:  Fk>at  100,  CM  IS,  23  November  1080,  0000  UTC 
3.0.1  Overview 

Tbe  drciiitoas  path  that  BAFOS  float  199  (Figure  29)  took  through  the  Central  Array 
wai  unique  for  SYNOP  floats.  The  thermocline  topography  (Figure  30)  shows  two  distinct 
features:  an  east-to-west  Gulf  Stream  to  the  north  of  a  ^clonic  cold  eddy.  Float  199  left 
the  stream  and  was  entrained  into  the  cold  eddy,  circling  three  times  before  being  expdled. 


Figure  29:  As  in  Figure  1,  but  for  float  199.  The  three  closest  RAFOS  fixes  to  Chf  15  (first 
approach)  ate  boxed.  Launch  date:  1258  UTC,  10  November  1989. 


3.9.2  Vertical  motion  results 

=  0.24  mm  s“*  ,  =  0.23  mm  s”*  ,  to,,,  =  -0.13  mm  s”* 

As  float  199  circled  the  cold  eddy  and  passed  close  to  15  (the  first  time),  it  was  moved 
from  588  to  593  to  574  db,  a  mean  upwdling  of  0.24  mm  s~^  (Figure  29);  at  this  time,  the 
float  was  in  the  middle  of  a  vertical  displacement  of  some  70  m.  CM  data  (Figure  31)  for 
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Figure  30:  As  in  Figure  2,  except  for  0000  UTC  23  November  1989.  The  track  of  float  199 
is  indicated  by  dots,  and  the  location  of  CM  15  is  marked  by  a  drde. 


the  same  time  also  suggests  upward  vertical  motion:  although  cnrreats  are  baddag  with 
decreadag  depth  (suggestive  of  downwdling)  (Figure  31a,c,e),  there  is  a  sfaaaltaaeoos  de¬ 
crease  ia  temperature  associated  with  a  rise  ia  isothermal  surfsoss  that  overwhelms  this 
desceat  (Figure  31b, d^)-  Vertical  motioa  as  diagaosed  from  lES  data  (Flgaie  30),  hoir- 
ever,  suggests  weak  dowawdliag  as  a  result  of  both  local  vortidty  teadoiqr  aad  vorUcity 
advectioa. 
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Fiiura  31:  As  ia  Figure  3,  except  for  CM  15  for  the  live  days  sarroending  float  IM’s  closest 
(first)  ^proadi  to  15,  which  time  is  indicated  by  vertical  line. 
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3.10  Case  10:  Float  109,  CM  IB,  28  Novambar  1089, 1600  UTC 
3.10.1  Ovarviaw 

Sea  Caaa  9.  At  this  time,  the  float  is  still  circling  arotind  the  cold  eddy  (Figure  32)  that 
flgufes  pimninently  ia  the  IBS  topography  (Figure  33). 


satismaumcB 


Figure  33:  As  ia  Figure  29,  but  for  sea»d  close  :4>proach. 


8.10.2  Vartkal  motloa  raeulta 

s  -0.13  mm  s"*  ,  »  0.04  mm  s~*  ,  -  0Ji4  mm  s~* 

As  RAFOS  float  199  passed  hf  CM  15  for  the  seamd  tioie,  it  aiofvad  firam  613  to  627 
to  621  db  (Figure  32),  a  meaa  dowaweHii^  of  -0.13  aua  s~*  ;  howevec,  hnawdlatily  after 
passi^  15,  the  float  upwdled  about  75  m  ia  2  dajrs.  CM  data  (Figure  34)  skomi  apmsBag, 
caused  mostly  by  local  tmaperature  chaagss  (Figure  34b,d4):  cumats  are  sudtddag  betwem 
bacldag  and  ve^ag  whea  the  float  passed  by,  so  wtical  motiem  from  tide  eftbet  was  aear 
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Figure  33:  A«  in  Figure  30,  except  for  1600  UTC  28  November  1969. 


zero  (Figure  34n,c,e).  lES  verticnl  motion  (Figure  33)  nt  the  znme  time  suggeeted  nperard 
motkm  (0.54  mm  )  forced  moetly  by  cyclonic  vortidty  advection. 
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3.11  Case  11:  Float  201,  CM  II,  23  November  1039,  0300  UTC 
3.11.1  Overview 

Float  201  moved  thoagh  the  aorthern  part  of  the  Central  Array  (Figure  35)  at  a  time  when 
the  thermocline  topography  (Figure  36)  indicated  a  vigorous  ccdd  ring  in  the  southern  part 
of  the  array  (In  fact,  float  199  was  circling  this  eddy  as  float  201  moved  through  the  array). 
zi2  data  also  show  the  sharp  trough  in  the  cmiter  of  the  Central  Array  whidi  float  201 
passed  just  before  being  expdled  from  the  stream  to  the  north  of  the  Cmitral  Array. 


Figure  35;  As  in  Figure  1,  but  for  float  201.  The  three  closest  RAFOS  fixes  to  CM  II  ate 
boxed.  Launch  date:  1550  UTC,  15  November  1989. 


3.11.2  Vertical  motion  results 

^KAr  “  »  ^CH  —  *”*  *  ®/*»  “ 

As  fioat  201  passed  by  CM  II,  it  was  just  be^nning  an  ascent  that  would  take  the  float 
firom  600  db  up  to  150  db  (Figure  35).  For  the  three  times  that  the  float  was  closest  to  II, 
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Figure  36:  As  in  Figure  2,  except  for  0800  UTC  23  November  1989.  The  trad(  of  float  201 
is  indicated  by  dots,  and  the  location  of  CM  II  is  marked  by  a  circle. 

it  moved  from  575  to  560  to  519  db,  a  mean  npwdling  of  0.97  mm  s"*  .  CM  data  for  this 
time  also  suggest  strong  upward  motion  (Figure  37).  Not  only  are  tempmatnres  cooling 
significantly  (suggestive  of  strong  upwelHng),  but  there  is  also  significant  veering.  Note  that 
this  is  a  case  where  the  computation  of  w^tt  a*ing  leads  to  an  (assumed)  erroneously 
large  component  of  vertical  motion  due  to  backing  and  veering  currents  (Figure  37a,c,e) 
that  does  not  occur  if  is  computed  using  v|*  <-  u^.  lES  topogr^^  also  suggests 
strong  upward  motion,  forced  entirely  by  the  cyclic  vortidty  advectkm  downstream  of 
the  sharp  trough  at  x  =  -40  (Figure  36b,d). 
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Figure  37:  As  in  Figure  3,  except  for  closest  approach  by  float  201  to  II,  which  time  is 
indicated  by  vertical  line. 
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3.12  Case  12:  Float  207,  CM  12,  4  September  1089,  0800  UTC 
S.12.1  Overview 

During  float  207’8  lengthy  passage  through  the  Central  Array  (Figure  38),  the  Central 
Array  was  characterized  by  unusually  intense  eddy  activity  (Figure  39).  As  a  result,  the 
float  followed  a  complex  route,  entering  the  Central  Array  four  times,  twice  from  the  eastern 
edge.  At  the  time  of  passing  12,  the  float  was  moving  to  the  west  on  the  southern  edge  of 
an  anticyclonic  (warm)  eddy. 


DAYS  ARB  UDNCS 
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Figure  38:  As  in  Figure  1,  but  for  float  207.  The  three  closest  RAFOS  Axes  to  CM  12  are 
boxed.  Launch  date:  1858  UTC,  23  August  1989. 


3.12.2  Vertical  motion  results 

=  -0.25  mm  s“^  ,  -  -0.06  mm  s“*  ,  w^gg  =  -0.24  mm 

Float  207  (Figure  38)  moved  from  158  to  164  to  172  db  during  the  time  it  was  closest  to 
CM  12,  at  a  level  in  the  ocean  that  was  the  shallowest  of  all  the  RAFOS  floats  considered  in 
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Figure  39:  As  in  Figure  2,  except  for  0800  UTC  4  September  1989.  The  track  of  float  207 
is  indicated  by  dots,  and  the  location  of  CM  12  is  marked  by  a  circle. 


this  report  (Float  201  moved  to  a  shallower  level  after  exiting  the  Central  Array,  however). 
Current  meter  data  at  the  same  time  also  suggested  downwelling  —  albeit  very  weak  — 
(Figure  40),  forced  by  backing  currents.  The  effect  of  the  backing  currents  dominates 
upwelling  diagnosed  as  a  result  of  cooling  temperatures.  lES  data  suggests  downwelling 
forced  mostly  by  local  vorticity  changes  (Figure  39c). 
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Figure  40;  As  in  Figure  3,  except  for  closest  approach  by  float  207  to  12,  which  time 
indicated  by  vertical  line. 


3.13  Case  13:  Float  209,  CM  15,  13  September  1989,  0800  UTC 
3.13.1  Overview 

RAFOS  float  209  moved  steadily  through  the  Central  Array  (Figure  41)  on  a  path  that 
suggests  the  float  moved  around  a  large  meander  trough  in  the  eastern  part  of  the  array. 
This  trough  is  obvious  as  well  in  the  zu  topography  (Figure  42),  which  also  shows  a  warm 
eddy  moving  through  the  slope  waters  to  the  north  of  the  stream. 
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Figure  41:  As  in  Figure  1,  but  for  float  209.  The  three  closest  RAFOS  Axes  to  CM  15  are 
boxed.  Launch  date:  1034  UTC,  7  September  1989. 


3.13.2  Vertical  motion  results 

=  -1.53  mm  s“^  ,  =  -2.14  mm  s~*  ,  =  0.22  mm  s“^ 

Float  209  (Figure  41)  moved  from  584  to  629  to  672  db  as  it  moved  closest  to  CM  15.  The 
mean  downwelling  from  this  motion  is  among  the  largest  observed  in  the  Central  Array.  At 
the  same  time,  CM  data  (Figure  43)  shows  currents  are  backing,  which  is  consistent  with 
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Figure  43:  As  in  Figure  3,  except  for  closest  approach  by  float  209  to  15,  which  time 
indicated  by  vertical  line. 


3.14  Case  14:  Float  210,  CM  H4,  6  October  1980,  0800  UTC 
3.14.1  Overview 

Float  210  passed  through  the  Central  Array  twice,  suggesting  significant  interaction  with 
nearby  rings  (Figure  44).  Indeed,  zu  data  (Figure  45)  suggests  that  fioat  210  is  circling  a 
developing  ring. 
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Figure  44;  As  in  Figure  1,  but  for  float  210.  The  three  closest  RAFOS  fixes  to  CM  H4  are 
boxed.  Launch  date:  1326  UTC,  26  September  1989. 


3.14.2  Vertical  motion  results 

vJka,’’  ~  “0*02  mm  s"*  ,  Wcm  =  -0.19  mm  s~*  ,  w,gg  —  -0.84  mm  s~^ 

As  float  210  moved  by  CM  H4  (Figure  44),  it  moved  only  from  636  to  644  to  637 
db,  despite  the  anticyclonic  curvature  to  the  path  that  might  suggest  strong  downweUing. 
Current  meter  data  at  the  same  time  also  reflects  the  slight  downward  motion  (Figure  46). 
Temperatures  are  warming,  which  suggests  downweUing  of  isotherms,  but  this  motion  is 
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Figure  45:  As  in  Figure  2,  except  for  0800  UTC  6  October  1989.  The  track  of  float  210  is 
indicated  by  dots,  and  the  location  of  CM  H4  is  marked  by  a  circle. 


mostly  offset  by  the  veering  current.  zi2  data  suggests  descent  forced  by  vorticity  advection 
(Figure  46). 
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Figure  46:  As  in  Figure  3,  except  for  closest  approach  by  float  210  to  H4,  which  time  is 
indicated  by  vertical  line. 
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3.15  Case  15:  Float  210,  CM  15,  21  October  1080,  0000  UTC 
3.15.1  Overview 

Float  210’s  second  passaj^e  through  the  central  array  occurred  as  the  float  (Figure  47), 
circulating  around  a  cold  eddy  (Figure  48),  entered  the  southern  edge  of  the  Cratral  Array 
and  then  turned  to  the  west. 
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Figure  47:  As  in  Figure  44.  The  three  closest  RAFOS  Axes  to  CM  15  are  boxed.  Launch 
date:  1326  UTC,  26  September  1989. 


3.15.2  Vertical  motion  results 

^RAF  ==  0-36  mm  s~*  ,  =  1.07  mm  s“*  ,  =  -0.33  mm  s“^ 

As  float  210  passed  by  CM  15,  the  float  upwelled  from  571  to  557  to  550  db,  (Figure  47) 
consistent  with  the  float’s  cyclonicaily  curved  path.  Indeed,  the  float  was  in  the  middle  ci  a 
three-day  upward  excursion  of  100  m  as  it  passed  15.  CM  data  also  suggest  strong  npwdling 
(Figure  49):  there  is  significant  veering  to  the  currents  and  isotherms  are  rising,  as  indicated 
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Figure  4^:  As  in  Figure  2,  except  for  0000  UTC  21  October  1989.  The  track  of  float  210  is 
indicated  by  dots,  and  the  location  of  CM  15  is  marked  by  a  circle. 


by  the  decreasing  temperatures  (Figure  49b).  lES  data,  however,  show  downweUing,  forced 
both  by  anticyclonic  vortidty  advection  and  by  local  tendency  (Figure  48).  We  note  again 
that  site  15  is  not  well  located  for  w,gg  computations. 
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Figure  49:  As  in  Figure  3,  except  for  closest  approach  by  float  210  to  15,  which  time  is 
indicated  by  vertical  line. 
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3.16  Case  16:  Float  211,  CM  M13,  8  October  1989,  0000  UTC 
3.16.1  Overview 

Float  211  (Figure  50)  moved  on  an  anticycloiiically  curved  path  the  entire  time  it  was  within 
the  Central  Array.  This  was  because  of  a  crest/trough  structure  in  the  rij  topography 
(Figure  51)  that  was  directing  the  main  Gulf  Stream  flow  out  the  southern  Central  Array 
border. 
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Figure  50:  As  in  Figure  1.  The  three  closest  RAFOS  Axes  to  CM  M13  are  boxed.  Launch 
date:  1243  UTC,  30  September  1989. 


3.16.2  Vertical  motion  results 

=  -0.44  mm  s“'  ,  =  -0.68  mm  s“*  ,  w,gg  =  -0.78  mm  s“* 

During  the  three  closest  positions  to  CM  M13,  RAFOS  float  211  was  at  pressure  levels 
613,  659,  and  638  db,  for  a  mean  downward  motion  of  -0.44  mm  s~^  (Figure  50).  CM 
data  (Figure  52)  also  suggested  downward  motion:  temperatures  at  M13  were  increasing. 
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Figure  51:  As  in  Figure  2,  except  for  0000  UTC  8  October  1989.  The  track  of  float  211  is 
indicated  by  dots,  and  the  location  of  CM  M13  is  marked  by  a  circle. 


reflecting  downward-moving  isotherms,  which  dominated  the  vertical  motion  diagnosed  by 
the  veering  currents  at  M13.  lES  data  also  suggests  downward  motion  (Figure  51),  as 
should  be  expected  given  the  strong  anticyclonic  vorticity  advection  over  the  site. 
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Figure  52:  As  in  Figure  3,  except  for  closest  approach  by  float  211  to  M13,  which  time  is 
indicated  by  vertical  Une. 
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3.17  Case  17;  Float  216,  CM  II,  26  November  1989,  0000  UTC 
3.17.1  Overview 

Float  216  (Figure  53)  moved  through  the  Central  Array  at  a  time  when  a  cold  eddy  was 
in  the  southeast  corner  of  the  array  (Figure  54);  however,  the  float’s  motion  through  the 
array  was  evidently  unaffected  by  this  cold  ring. 
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Figure  53:  As  in  Figure  1.  The  three  closest  RAFOS  fixes  to  CM  11  are  boxed.  Launch 
date:  1034  UTC,  18  November  1989. 


3.17.2  Vertical  motion  results 

=  -0.02  mm  s“^  ,  Wcm  =  0.14  mm  8~*  ,  w,gs  =  0.86  mm  s“* 

Float  216  passed  by  CM  II  just  before  downwelling  more  than  200  m  in  a  week  (Fig¬ 
ure  53).  Close  to  II,  however,  the  float  moved  from  567  to  576  to  568  db  on  a  path  the 
trajectory  of  which  was  strmght.  lES  data  (Figure  54)  shows  significant  curvature  to  the  zia 
field,  however,  with  cyclonic  vorticity  advection  forcing  upward  motion  at  II  (in  fact,  the 
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Figure  54:  As  in  Figure  2,  except  for  0000  UTC  26  November  1989.  The  track  of  float  216 
is  indicated  by  dots,  and  the  location  of  CM  II  is  marked  by  a  circle. 


fields  look  very  similar  to  those  associated  with  float  201).  CM  data  for  the  same  time  also 
suggest  a  Gulf  Stream  that  is  not  straight  (there  is  time  variability  in  the  direction  (Fig¬ 
ure  55a)).  Upward  motion  is  forced  by  veering  currents,  overwhelming  the  weak  downward 
motion  associated  with  the  isotherm’s  descent  (Figure  55b,d). 
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3.18  Case  18:  Float  221,  CM  12,  6  January  1990,  0800  UTC 
3.18.1  Overview 

Float  221  (Figure  56)  moved  quickly  through  the  Central  Array  on  a  straight  path, 
data  (Figure  57)  also  shows  a  mostly  straight  Gulf  Stream  path  (in  which  are  embedded 
minor  troughs/crests)  with  a  cold  eddy  far  to  the  south  of  the  stream. 


Longitude  (W) 

Figure  56:  As  in  Figure  1.  The  three  closest  RAFOS  fixes  to  CM  12  are  boxed.  Launch 
date:  0253  UTC,  30  December  1989. 


3.18.2  Vertical  motion  results 

^KAF  —  »  ^/ES  ~ 

Float  221  moved  from  490  to  483  to  480  db  in  the  time  surrounding  its  closest  passage  to 
CM  12  (Figure  56).  lES  data  suggest  this  upward  ascent  is  forced  both  by  cyclonic  vorticity 
advection  and  by  local  vorticity  tendencies  (Figure  57).  CM  data  (Figure  58)  also  reflect 
the  upward  motions,  as  there  is  significant  veering  motion  and  slight  cooling. 
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Figure  57:  As  in  Figure  2,  except  for  0800  UTC  6  January  1990.  The  track  of  float  221  is 
indicated  by  dots,  and  the  location  of  CM  12  is  marked  by  a  circle. 
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3.19  Case  10:  Float  224,  CM  H4,  20  January  1090,  0800  UTC 
3.19.1  Overview 

Float  224  (Figure  59)  passes  through  the  Central  Array  at  a  time  when  a  cold  eddy  is 
apparently  interacting  with  the  Gulf  Stream  (Figure  60).  This  interaction,  however,  does 
not  seem  to  affect  the  float. 


Longitude  (W) 

Figure  59:  As  in  Figure  1.  The  three  closest  RAFOS  fixes  to  CM  H4  are  boxed.  Launch 
date:  1200  UTC,  15  January  1990. 


3.10.2  Vertical  motion  results 

=  -0.18  mm  s“*  =  -0.10  mm  s“’  =  -0.30  mm  s“^ 

As  it  passed  close  to  CM  H4,  float  224  moved  from  477  to  497  to  487  db.  The  pressure 
trace  (Figure  59)  suggests  that  the  float  downwelled  until  0800  UTC  20  January  1990  (when 
it  was  closest  to  H4),  and  then  it  leveled  off,  or  upwelled  slightly.  lES  data  suggests  this 
may  also  be  true;  although  downwelling  is  diagnosed  in  the  vicinity  of  H4,  upwelling  is 
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Figure  60:  As  in  Figure  2,  except  for  0800  UTC  20  January  1990.  The  track  of  float  224  is 
indicated  by  dots,  and  the  location  of  CM  H4  is  marked  by  a  circle. 


occurring  just  downstream,  both  forced  primarily  by  vorticity  advection  (Figure  60).  CM 
data  also  sugge'^t  downwelling;  veering  currents  are  forcing  upward  motion  (Figure  61a,c,e) 
that  is  offset  by  the  greater  downward  motion  associated  with  local  warming/descending 
isotherm  (Figure  61b,d,f).  This  is  another  example  where  the  temperature  tendency  at  700 
m  is  much  stronger  than  the  average  of  the  400  m  and  1000  m  tendendes.  The  1000  m  and 
400  m  sensors  are  respectively  below  and  above  the  main  thermocline. 
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4  Comparing  the  three  different  methods 

4.1  Statistics 

What  are  the  general  statistics  of  vertical  motion  in  the  Central  Array  from  the  three  data 
sources?  Table  1  shows  the  means,  standard  deviations,  maxima,  minima,  and  the  number 
of  observations  for  the  two-year  period  for  each  data  source.  All  data  sources  indicate  mean 
downweOing,  which  is  consistent  with  expectations  from  quasi-geostrophic  motion  and  the 
time  average  of  equation  (7)  given  the  the  observed  crest  just  upstream  of  the  Central  Array 
(near  the  Gulf  Stream  *node’  at  70*  W)  and  a  trough  near  the  eastern  boundary  of  the  array 
in  the  mean.  The  larger  mean  value  of  may  reflect  a  sampling  bias:  floats  remained 
mostly  within  the  strongly  bsuoclinic  structure  of  the  central  Gulf  Stream  when  traversing 
the  Central  Array,  whereas  CM  and  lES  measured  vertical  motions  throughout  the  r^on 
spanning  also  the  Slope  and  Sargasso  waters  where  large-scale  vertical  motions  are  smaller. 
Note  that  downwelling  averaged  over  all  tracks,  including  those  that  escaped  the 
Gulf  Stream,  is  an  order  of  magnitude  smaUer  than  within  the  Central  Array.  We  are 
examining  a  segment  of  the  Gulf  Stream  that,  for  these  26  months  at  least,  experienced 
mean  downwelling. 

Ihble  1:  Vertical  Motion  Statistics:  Mean,  maximum,  and  minimum  vertical  motion  (in 
mm  s"*)  with  standard  deviation  and  number  of  (not  all  independent)  observations  for 
CM,  lES,  RAFOS  floats  in  the  Central  Array  only,  and  for  all  RAFOS  floats 


Source 

Mean 

Std.  Deviation 

Maximum 

Minimum 

#obe 

Current  Meters 

-0.015 

0.347 

3.676 

-2.297 

25419 

Mapped  lES 

0.440 

-3.794 

188640 

RAFOS  in  central  array 

0.617 

2.10 

-2.66 

947 

AU  RAFOS 

0.532 

3.30 

-3.34 

5939 

To  test  the  statistical  ugniflcance  of  mean  downwelling  in  Table  1,  we  computed  the 
uncertainty  of  the  mean  for  timeseries  of  at  each  CM  site,  and  at  the  ten  extrema  in 
the  mean  w,gg  Add  in  Figure  62  using  techniques  in  Bendat  and  Piersd  (1973)  and  Dewar 
and  Bane  (1985).  For  each  of  the  ten  extrema,  the  mean  value  differed  from  zero  by 
at  least  one  standard  deviation  (and  usually  two).  Mean  values  of  tOcg  also  wme  usually 
at  least  one  standard  deviation  from  zero,  except  for  several  sites  where  Wcu  was  very 
nearly  zero.  The  mean  downwelling  is  associated  with  the  fact  that  the  lES  instruments 
preferentially  sampled  a  crest-to-trongh  segment  of  the  Gulf  Stream.  It  follows  from  the 
work  of  Bower  (1989)  that  downwelling  should  occur.  How  representative  of  the  long-term 
mean  are  our  results?  Lee  (1993)  shows  a  trough  in  the  Central  Array  for  6  out  8  years 
using  infrared  satellite  data.  We  do  believe  the  values  for  mean  vertical  motion  in  Table  I 
have  some  statistical  significance. 

4.2  Case  Studies 

We  now  present  several  Figures  that  give  evidence  that  the  vertical  motion  determined  from 
ail  three  data  sources  is  consistent.  Figure  63  shows  two  time-series  of  Wcm  and  w,gg  at 
two  unrelated  sites  and  times,  G2  and  II.  There  are  clear  similarities  between  and 
v>,gg:  both  show  events  of  up-  and  downwelling  that  develop  and  pass  coherently  with  some 
high-frequency  noise,  extrema  are  approximately  colocated,  signal  strength  is  similar  in  the 
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Figure  62:  Mean  zn  topography  (bold  lines,  contoured  every  100  m)  and  mean  vertical 
motion  from  (7)  (thin  lines  contoured  every  0.05  mm  negative  values  dashed)  for  the 
period  15  June  1988  -  7  August  1990  within  the  Central  Array.  P(^t  values  (±  standard 
deviation)  of  extrema  as  indicated,  x-  and  y-  axis  labels  are  in  kilometers. 


two  lines,  and  the  episodic  nature  is  obvious.  We  have  made  a  similar  comparistm  for  all  CM 
sites  in  the  two-year  deployment  period;  the  similarity  evident  in  Figure  63  characterises  all 
sites:  large  vertical  motion  in  one  data  source  usually  has  a  corresponding  extremum  in  the 
other;  however,  some  counterexamples  will  be  shown.  Elxamples  of  the  w  records  along  the 
line  of  current  meters  (G2  and  G3  for  both  Year  1  and  Year  2),  along  the  'H’  line  (H2, 
H3,  H4,  and  H6  for  both  Year  1  and  Year  2,  H5  for  only  Year  1),  along  the  ‘F  line  (II,  12, 
13, 14, 15  for  both  Year  1  and  Year  2),  and  at  M13  (Year  2)  are  shown  in  Figures  64 —  76. 
In  each  of  these  time  series,  the  characteristics  highlighted  in  Figure  63  are  evident.  In  each 
time  series,  a  large  peak  in  vertical  motion  in  one  method  generally  has  a  corresptmding 
peak  in  the  other  method,  although  some  time  series  (e.y..  Figure  69)  are  better  than  others 
(Figure  75).  There  is  inter-time  series  variability  as  well.  Compare,  for  instance.  Year  1 
in  Figure  68,  which  shows  little  coherence,  with  Year  2,  which  shows  considerably  more. 
Despite  occasional  short-period  discrepancies  between  methods,  we  show  below  that  time 
series  of  and  w,gg  are  coherent  at  periods  longer  than  12  days. 
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Figure  63:  Time  series  of  vertical  motion.  (A)  CM  G2  and  (B)  CM  II  witk  dates  as  indi¬ 
cated  on  the  z-axis.  Vertical  motion  from  CM  data  using  equation  (4)  in  indicated 
with  the  dotted  Une;  that  from  lES  data  {w,gs)  &nd  equation  (7)  is  indicated  with  the  solid 
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Figure  64:  Time  series  of  w,gg  (solid  line)  and  Wcu  (dashed  line)  for  CM  G2  for  Year  1 
and  Year  2. 
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Figure  65:  As  in  Figure  64,  but  for  G3. 
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Figure  67:  As  in  Figure  64,  but  for  H3. 
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Figure  68:  As  in  Figure  64,  but  for  H4. 
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Figure  70:  As  in  Figure  64,  but  for  H6. 
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Figure  71:  As  in  Figure  64,  but  for  II. 
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Table  2  presents  comparisons  of  vertical  motion  for  all  23  cases  in  which  a  float  passed 
within  10  km  of  a  CM  mooring  with  working  instrumentation  at  400,  700,  and  1000  m. 
These  are  the  cases  discussed  in  Section  3.  Figure  77a  is  a  scatterplot  of  RAFOS  vertical 
motions  versus  CM  vertical  motions,  with  data  taken  from  Table  2.  A  linear  relationship 
is  clearly  evident  with  correlation  =  0.82.  That  the  slope  of  the  line  is  greater  than 
unity  could  arise  from  our  choices  of  values  for  or  or  from  being  biased  low. 
However,  the  plot  underscores  the  excellent  agreement  between  and 


Table  2:  Vertical  Motion  Comparisons:  Vertical  motions  from  RAFOS  float,  CM,  and  lES 
data  Wca/»  “'/es'  respectively),  with  units  of  mm  s~'  for  date/time  shown  and 

for  RAFOS  float  and  CM  indicated.  Depth  of  the  RAFOS  float  at  the  comparison  time  (in 
meters)  is  also  shown. 


Float 

CM 

date/time 

EHfl 

I!3i9l 

Depth 

RAF123 

H4 

4Jul88/00z 

-0.58 

-0.98 

-1.42 

550.3  m 

RAF129 

14 

23Sep88/08z 

0.19 

0.53 

0.04 

710.7  m 

RAF129 

14 

23Set.C8/16z 

0.72 

0.48 

0.05 

706.4  m 

RAF136 

14 

8Dec88/16z 

0.81 

0.43 

0.17 

668.4  m 

RAF141 

14 

2Nov88/16z 

-0.87 

-0.94 

-0.43 

603.0  m 

RAF175 

14 

20Jan89/16z 

-0.56 

-0.05 

-0.63 

363.7  m 

RAF176 

15 

7Feb89/00z 

-0.04 

0.02 

0.25 

761.5  m 

RAF176 

15 

7Feb89/08z 

-0.05 

0.03 

0.26 

763.0  m 

RAF176 

15 

7Feb89/16z 

0.04 

0.02 

0.27 

764.5  m 

RAF176 

15 

8Feb89/00z 

-0.03 

-0.05 

0.26 

760.7  m 

RAF  194 

15 

15Nov89/00z 

0.71 

0.23 

Bim 

662.2  m 

RAF194 

14 

16Nov89/00z 

0.86 

1.20 

-1.06 

584.3  m 

RAF  199 

15 

23Nov89/00z 

0.24 

0.23 

593.2  m 

RAF199 

15 

28Nov89/16z 

-0.13 

0.04 

0.54 

627.4  m 

RAF201 

11 

23Nov89/08z 

0.97 

1.51 

1.66 

560.4  m 

RAF207 

12 

4Sep89/08z 

-0.25 

-0.24 

164.0  m 

RAF209 

15 

13Sep89/08z 

-1.53 

-2.14 

0.22 

629.2  m 

RAF210 

H4 

60ct89^8z 

-0.02 

-0.19 

642.4  m 

RAF210 

15 

210ct89/00z 

0.36 

1.07 

-0.33 

557.4  m 

RAF211 

M13 

80ct89/00z 

-0.44 

-0.68 

-0.78 

658.9  m 

RAF216 

11 

26Nov89/00z 

-0.02 

0.14 

0.86 

576.2  m 

12 

6Jan90/08z 

0.18 

0.21 

0.45 

482.9  m 

RAF224 

H4 

20Jan90/08z 

-0.18 

-0.10 

-0.30 

496.7  m 

The  comparison  between  and  w,ps  for  the  same  23  cases  (Figure  77b)  shows 

more  scatter,  in  part  because  the  cases  were  selected  for  floats  passing  CM  sites;  several 
of  these  cases  are  not  in  the  best  lES  mapping  region.  Nevertheless,  after  excluding  the 
two  obvious  outliers  that  are  discussed  below,  the  correlation  =  0.32,  which  is  significant 
at  a  confidence  level  of  greater  than  99%.  That  the  slope  of  the  line  is  less  than  unity 
suggests  there  may  be  a  bias  towards  underestimating  w,pg  that  would  result  from  the 
spatial  filtering  of  the  lES  maps. 

Careful  scrutiny  of  Table  2  shows  that  Wj^g  consistently  is  of  the  wrong  sign  at  15. 
The  computation  of  vorticity  and  vorticity  advection  at  site  15  is  adversely  impacted  by 
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that  site’s  close  proximity  to  edge  of  the  domain.  After  excluding  all  from  site  15  in 
Table  2,  we  plot  the  remaining  values  in  Figure  77c,  from  which  we  compute  a  correlation 
of  0.51,  significant  at  a  confidence  level  exceeding  90%.  To  understand  the  causes  of  errors 
in  Figures  77b,c,  we  have  closely  investigated  the  vertical  motion  fields  associated  with  the 
one  remaining  outlier  in  Figure  77c.  Two  possible  error  sources  are  discussed  below. 


-3.0-1.S  0.0 
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Figure  77:  (A)  Scatterplot  of  versus  tv^^,  for  data  as  listed  in  Table  2.  (B)  Scatterplot 
of  versus  as  listed  in  Tabic  2.  (C)  As  in  (B),  but  excluding  W/gg  at  site  15.  The 
box  point  in  (B)  and  (C)  is  discussed  further  in  text. 


The  point  in  the  lower  right  quadrant  of  Figure  77b,c  is  from  the  case  when  BAFOS 
float  194  passed  by  CM  14  at  0000  UTC  16  November  1989  (See  Section  3.8).  Both  CM 
and  float  vertical  velocities  are  upward;  however,  w,gg  is  strongly  downward.  This  may  be 
caused  by  poor  lES  data  at  site  I4  during  winter  of  Year  2,  when  this  lES  suffered  from 
many  bad  data  returns,  necessitating  frequent  interpolations  to  compute  thermocline  depth 
at  that  site  (Fields  and  Watts  1991).  Such  data  degradation  would  of  course  have  a  negative 
impact  on  the  computed  streamfunction  field  and  all  quasigeostrophic  calculations  near  14 
for  winter  of  1989-1990  data. 

It  is  far  more  likely  that  errors  were  introduced  into  w,gg  at  this  time  because  significant 
currents  existed  at  deep  levels.  lES  measurements  determine  the  baroclinic  component  of 
the  flow  only.  IF  there  is  a  significant  barotropic  flow,  then  the  vorticities  and  advections 
computed  in  (7)  using  only  lES  data  will  not  be  representative  of  the  true  flow.  Figure  27 
presents  the  lES  topography  and  w,gg  for  this  time,  as  well  as  the  track  of  float  209  as  it 
traversed  the  central  array.  The  deep  currents  at  the  CM  sites  are  also  plotted.  Note  the 
unusually  strong  flow  at  3500  m  at  14.  In  this  area,  strong  upward  motion  associated  with 
the  cyclonically  curved  onshore  barotropic  flow  overwhelms  the  downward  motion  associated 
with  the  baroclinic  flow.  In  fact,  this  event  is  associated  with  the  strongest  deep  currents 
for  the  Wjgg  cases  plotted  in  Figure  77c:  |t^5ool  =  0.19m  s“*,  which  is  a  value  nearly  two 
standard  deviations  above  the  mean  speed  at  14  during  Year  2.  Note  that  strong  deep 
currents  do  not  necessarily  mean  poor  values  of  u),gg:  Case  11  (Section  3.11)  also  had  deep 
velocities  in  excess  of  0.15  m  s~^  and  there  is  good  agreement  between  and  tB,gg.  For 
this  case,  the  deep  currents  flow  parallel  to  the  Gulf  Stream;  for  the  case  above,  the  deep 
flow  crosses  under  the  Gulf  Stream.  Where  the  barotropic  flow  is  significant,  especially  if  it 
is  spatially  non-uniform,  errors  in  w,gg  are  quite  likely.  Because  it  is  possible  to  account  for 
the  deep  currents,  however,  (for  example  by  combining  a  deep  barotropic  streamfunction  — 
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measured  by  current  meters  or  perhaps  by  electric  field  measurements  —  with  the  baroclinic 
streamfunction)  we  do  not  consider  this  an  unbearable  problem.  Indeed,  this  is  the  subject 
of  an  ongoing  followup  study. 

4.3  Coherence  between  and 

The  time-series  estimates  of  and  to,g;g  (e.g..  Figures  63  and  64-  76)  show  smoothly 
varying  signals  with  varying  degrees  of  high-frequency  noise  superimposed.  We  spectrally 
decomposed  each  pair  of  records  and  computed  their  statistical  coherence  at  all  sites  along 
the  central  mooring  H-line  (whose  time  series  are  shown  in  Figures  66-  70)  at  those  CM 
sites  that  had  all  three  upper  CMs  functioning  and  that  were  surrounded  by  functioning 
lESs  (these  were  H3,  H5,  and  H6  from  Year  1  and  H4  and  H6  from  Year  2).  These  coherence 
results  were  averaged  for  the  whole  line  and  are  shown  in  Figure  78.  The  average  coherence 
between  the  two  independent  estimates  of  vertical  motion  is  greater  than  0.5  (at  far  above 
the  99%  confidence  level)  throughout  all  periods  longer  than  16  days.  At  periods  shorter 
than  12  days,  the  coherence  is  insignificant.  This  confirms  the  visual  impression  from 
Figures  63-  76  that  the  large  amplitude,  low-frequency  variations  of  w  from  both  CM  and 
lES  techniques  are  very  consistent,  while  the  snialler-amplitude,  high-frequency  variations 
are  noise-dominated.  In  Figures  79-  91,  we  show  the  coherence  for  all  of  the  sites,  calculated 
separately  for  Year  1  and  Year  2  because  the  instrument  arrays  changed  between  the  two 
years.  In  particular,  the  mapping  ability  of  the  lES  array  to  estimate  terms  in  the  vorticity 
equation  could  have  changed  between  the  two  years  -  as  it  appears  to  have  done  at  some 
sites  like  G2  (Figure  79)  and  H2  (Figure  81). 


Period  (days) 

Figure  78:  Coherence  as  a  function  of  frequency  for  and  for  selected  sites  on  the 
‘H’  line. 


Given  these  statistical  coherences,  how  well  can  we  use  equations  (4)  and  (7)  to  infer 
vertical  motion  using  data  from  the  CM  and  lES  arrays  that  were  designed  for  mesoecale 
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Figure  79:  Coherence  as  a  function  of  frequency  for  v),gs  and  for  CM  G2  Year  1  (solid 
line)  and  Year  2  (dotted  line). 


Figure  80:  As  in  Figure  79,  but  for  G3. 


sampling?  Both  techniques  involve  temporal  and  spatial  differentiation  of  measured  quan¬ 
tities,  an  analysis  which  accentuates  measurement  and  sampling  errors.  The  mors  in 
from  (4)  are  not  dominated  by  measurement  error  but  rather  by  submesoscale  features  for 
which  the  signal  is  incoherent  between  our  moored  u,  v,  T  sensors.  Errors  in  from 
(7)  are  dominated  by  instrument/analysis  error  in  the  streamfunction  Adds.  The  objective 
mapping  and  measurement  technique  filters  out  submesoscale  features,  induding  some  that 
may  be  sampled  at  individual  CM  moorings.  Neverthdess,  with  dther  technique,  the  major 
mesoscale  features  in  w  are  dear.  The  coherence  between  tn,,,  and  Wca  shown  in  Figure 
10  and  the  particularly  good  agreement  between  and  illustrated  in  Table  2  and 
Figure  77a,  t.e.  agreements  between  three  independent  techniques  for  measuring  w,  show 
that  the  V/c5(^>y«0  Adds  do  characterize  mesoscale  vertical  motions  in  the 

Gulf  Stream. 
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Figure  81:  As  in  Figure  79,  but  for  H2. 
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Figure  82:  As  in  Figure  79,  but  for  H3. 
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Figure  83:  As  in  Figure  79,  but  for  H4. 


Figure  84:  As  in  Figure  79,  but  for  H5  {Year  1  only). 
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Figure  85:  As  in  Figure  79,  but  for  H6. 
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Figure  86:  As  in  Figure  79,  but  for  II. 

1 


Figure  87:  As  in  Figure  79,  but  for  12. 


Figure  88:  As  in  Figure  79,  but  for  13. 


Figure  89:  As  in  Figure  79,  but  for  I4. 
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Figure  90:  As  in  Figure  79,  but  for  15. 
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Figure  91:  As  in  Figure  79,  but  for  M13. 
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5  Summary 

We  have  computed  and  compared  vertical  motions  from  three  different  sources  in  the  Gulf 
Stream  downstream  of  Cape  Hatteras.  Magnitudes  of  vertical  motion  as  measured  by 
RAFOS  floats  amd  as  estimated  from  CM  data  and  lES  data  are  1-2  mm  s~^,  with  rare 
values  of  up  to  3  mm  s~*.  Both  the  CM  and  lES  vertical  motions  are  smoothly  varying 
fields  not  dominated  by  the  submesoscale  and  measurement  noise  that  is  present.  It  is  easy 
to  track  features  from  day  to  day  in  both  datasets;  evolution  of  the  fields  occurs  only  slowly 
as  mesoscale  jet  features  evolve. 

CM  vertical  motions  agree  closely  with  observations  from  RAFOS  floats.  There  is  also 
considerable  agreement  between  lES  and  RAFOS  float  or  CM  vertical  motions,  and  very 
good  coherence  at  long  timescales  between  and  timeseries.  Indeed,  for  timescales 
of  greater  than  16d,  and  w,gg  show  coherence  >  0.5  at  far  above  the  99%  confidence 
level,  but  such  coherence  is  absent  for  timescales  of  less  than  12d.  This  lack  of  coherence  may 
be  caused  in  part  by  errors  related  to  the  presence  of  strong  deep  (barotropic)  velocities 
under  the  meandering  Gulf  Stream.  By  combining  the  lES  data  with  an  independent 
measure  of  the  barotropic  velocity,  it  should  be  possible  to  account  for  the  vertical  vdodty 
driven  by  deep  currents. 

To  our  knowledge,  this  is  the  first  study  in  which  vertical  motions  associated  with 
mesoscale  processes  in  the  ocean  have  been  observed  and  verified  by  independent  measure¬ 
ments  and  independent  dynamical  methods.  The  spatial  and  temporal  structures  of  the  w 
field  can  be  diagnosed  with  both  CM  and  lES  arrays;  the  Gulf  Stream  clearly  exhibits  sec¬ 
ondary  circulations  on  the  mesoscale  that  are  consistent  with  quasi-geostrophic  dynamics. 
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